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REMARKS 

Applicants respectfully request reconsideration of the present application in view of 
the foregoing amendments and in view of the reasons that follow. 

Claims 1,5, 12 and 16-18 are requested to be cancelled without prejudice or 
disclaimer. 

Claims 2-4, 6, 9-1 1 and 13-15 are currently being amended. 
Claims 19-21 are being added. 

This amendment adds, changes and deletes claims in this application. A detailed 
listing of all claims that are, or were, in the application, irrespective of whether the claim(s) 
remain under examination in the application, is presented, with an appropriate defined status 
identifier. 

After amending the claims as set forth above, claims 2-4, 6-11 and 13-21 are now 
pending in this application. 

Drawings 

The drawings were objected to. Specifically, the Office Action indicated that the low 
pass filter must be shown or the feature canceled from the claims. Applicants respectfully 
traverse this objection. The low pass filter in the claims is a mathematical function, not a 
structure, and as such need not be illustrated in the drawings. Moreover, at least Figure 5 
illustrates "low pass filtering" in step S50. 

Disclosure 

Paragraph 4 on page 3 of the Office Action stated "Applicant is reminded to use 
consistent language throughout the disclosure in order to facilitate finding support for the 
recited limitations, as well as to provide proper antecedence for all claim limitations," and 
listed examples a) through d). With respect to a), the claims have been amended to eliminate 
the acronym "SOC" after "charge rate." Applicants note that the present specification equates 
the charge rate with SOC. The claims have been amended to address the issue raised with 
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respect to b). With respect to c), applicants submit that the recitation in claim 3 of 1/Gl(s) as 
a "transfer function" is consistent with referring to 1/Gl(s) as being a low pass fiher or having 
low pass filter characteristics. With respect to d), applicants note that upper case "K" is used 
for intemal resistance while lower case "k" is used for time. Thus, different symbols are used 
to represent intemal resistance and time, and no confusion should result therefi-om. 

Claim objections 

The claims were objected to for informalities. The claims have been amended to 
address the issues raised in the Office Action, and applicants submit that the objections have 
been overcome. 

Allowable subject matter 

Applicants appreciate the indication that claims 4-11 and 15 contain allowable subject 

matter. 

Rejections under 35 U.S.C. § 103 

Claims 1-3, 12-14, and 16-18 were rejected under 35 U.S.C. § 103(a) as being 
unpatentable over "RESEARCH AND DEVELOPMENT OF THE MODEL-BASED 
BATTERY STATE OF CHARGE INDICATOR" to Torikai et al (Torikai). Claims 16-18 
were rejected under 35 U.S.C. § 103(a) as being unpatentable over "MODEL-BASED 
PREDICTIVE DIAGNOSIS FOR ELECTROCHEMICAL ENERGY SOURCES" to 
Kozlowski et al (Kozlowski). These rejections are moot with respect to claims 1,12 and 16- 
1 8, which have been canceled. Insofar as these rejections can be applied to the claims as 
gunended, applicants respectfully traverse for at least the following reasons. 

Claim 2 has been amended to be in independent form. In claim 2, equation (2) is 
substituted into battery model equation (1) to obtain equation (4), which does not include any 
offset term Vo . Thus, equation (1) can be converted into a product-and-addition equation, 
equation (4), between a measurable terminal voltage V, current I, and a parameter. This 
product-and-addition equation provides advantages in processing. Specifically, this allows an 
adaptive digital filter equation to be directly applied in a continuous time series. As a result, 
parameters including parameter h can be estimated in a batch processing manner. Open- 
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circuit voltage Vo may then be obtained using h among the batch-estimated parameters. 
Therefore, the charge rate is accurately estimated from "open-circuit voltage Vo" and a 
predetermined relationship between the open-circuit voltage and the charge rate of the 
secondary battery. 

In a similar fashion to independent claim 2, in independent claim 3, equation (2) is 
substituted into battery model equation (1) to obtain equation (4), which does not include an 
offset term Vq. Thus equation (1) can be converted into the product-and-addition equation, 
equation (4), between a measurable terminal vohage V, current I and a parameter. In a 
similar fashion to claim 2, in claim 3, this product-and-addition equation provides advantages 
in processing. Specifically, this allows an adaptive digital filter equation to be directly 
applied in a continuous time series. As a result, parameters A(s), B(s) and C(s) can be 
estimated, which are substituted into equation (5) of claim 3, which does not include an 
integration term of estimated parameters. Open-circuit voltage Vo may thus be obtained 
without effect of integration error. Therefore, the charge rate is accurately estimated from 
"open-circuit voltage Vq" and a predetermined relationship between the open-circuit voltage 
and the charge rate of the secondary battery. 

Torikai discloses that SOC is estimated by battery model equation (1) expressed with 
terminal voltage V, current I and SOC, 0. In contrast to claim 2, however, the battery model 
equation is not a product-and-addition equation between a measurable terminal voltage V, 
current I, and a parameter, as equation (2) of Torikai shows. As can be seen, the right hand 
side of equation (2) of Torikai is not a product-and-addition equation between current I and a 
parameter. Thus, equation (2) cannot be easily applied to an adaptive digital filter formula, in 
contrast to the equations of claims 2 or 3. 

Claims 13 and 14 are method claims corresponding respectively to claims 2 and 3, and 
are patentable for analogous reasons. 

The dependent claims are patentable for at least the same reasons as their respective 
independent claims, as well as for further patentable features recited therein. 
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New dependent claims 20 and 21 have been added, and are likewise believed to be 
patentable. 

The dependent claims are patentable for at least the same reasons as their respective 
independent claims, as well as for further patentable features recited therein. 

Applicants beUeve that the present application is now in condition for allowance. 
Favorable reconsideration of the application as amended is respectfully requested. 

The Examiner is invited to contact the undersigned by telephone if it is felt that a 
telephone interview would advance the prosecution of the present application. 

The Commissioner is hereby authorized to charge any additional fees which may be 
required regarding this application under 37 C.F.R. §§ 1.16-1.17, or credit any overpayment, 
to Deposit Account No. 1 9-074 L Should no proper payment be enclosed herewith, as by a 
check being in the wrong amount, unsigned, post-dated, otherwise improper or informal or 
even entirely missing, the Commissioner is authorized to charge the unpaid amount to 
Deposit Accoimt No. 19-0741. If any extensions of time are needed for timely acceptance of 
papers submitted herewith. Applicant hereby petitions for such extension under 37 C.F.R. 
§1 .136 and authorizes payment of any such extensions fees to Deposit Account No. 19-0741. 



Date ^7^^\k^?L.cM ^ > ><^C:> 

FOLEY & LARDNER LLP 
Customer Number: 22428 
Telephone : (202) 672-54 1 4 
Facsimile: (202) 672-5399 



Respectfully submitted. 




Richard L. Schwaab 
Attomey for Applicant 
Registration No. 25,479 

Thomas G. Bilodeau 
Attomey for Applicant 
Registration No. 43,438 
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Atty. Dkt. No. 023971-0333 

APPARATUS AND METHOD FOR ESTIMATING CHARGE RATE OF 

SECONDARY CELL 

BACKGROUND OF THE INVENTION: 

5 

Field of the Inven'tion 

[0001] The present invention relates to apparatus 

and method for estimating a charge rate (abbreviated 
10 as SOC) of a secondary cell. 

Descrip'tlon of the related art 

[0002] Japanese Patent Application First 

15 Publications No. 2000-323183 published on November 24, 
2000 No. 2000-268886 published on September 29, 2000, 
and a Japanese Paper titled Estimation of Open 
Voltage and Residual Values for Pb Battery by 
Adaptive Digital Filter announced by a Japanese 

20 Electrical Engineering Society (T.IEEE Japan), Volume 
112-C, No. 4, published on 1992 exemplify previously 
proposed SOC estimating apparatus for the secondary 
cell. That is to say, since the charge rate (or 
called State Of Charge, i.e., SOC) of the secondary 

25 cell has a correlation to an open-circuit voltage Vq 
(cell terminal voltage when its power supply of the 
cell is turned off, also called electromotive force 
or open voltage), the charge rate can be estimated 
when open voltage Vq is obtained. However, a 

30 considerable time is needed until the terminal 
voltage is stabilized after the power supply is 
turned off ( charge-and-di schar ge is ended) . Hence, a 
predetermined time duration is needed from a time at 
which the charge-and-discharge is ended to determine 
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an accurate open-circuit voltage Vq . Therefore, since 
immediately after or during the char ge /di s charge time 
or charge-and-discharge , it is impossible to 
determine an accurate open-circuit voltage and the 
5 charge rate cannot be obtained using the above- 
described method. Nevertheless, to determine the 
open-circuit voltage Vo, the open-circuit voltage Vq 
is estimated using a method disclosed in the above- 
described Japanese Patent Application First 
10 Publication No. 2000-323183. 

SUMMARY OF THE INVENTION: 
[0003] However, in the above-described method 

disclosed in the Japanese Patent Application 

15 Publication No. 2000-323183, open-circuit voltage Vq 
is calculated from a non-recursive (non-regression 
type) cell model (a model whose output value is 
determined only from a present value and past value 
of an input value) whose characteristic is wholly 

20 different from a physical characteristic of the cell 
for which an adaptive digital filter (sequential type 
model parameter identification algorithm) is used. 
The charge rate ±SOC)_ is used from this value. Hence, 
in a case where this method is applied to the actual 

25 cell characteristic (input: current, output: voltage), 
according to the cell characteristic, an estimation 
calculation is wholly converged or does not converge 
to a real value. Hence, it is difficult to estimate 
the charge rate (SQC) accurately. 

30 [0004] It is, hence, an object of the present 

invention to provide apparatus and method for 
estimating accurately the charge rate (SOC) for the 
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secondary cell and accurately estimating other 
parameters related to the charge rate (SOC) • 
[0005] According to one aspect of the present 

invention, there is provided a charge rate estimating 
5 apparatus for a secondary cell, comprising: a current 
detecting section capable of measuring a current 
flowing through the secondary cell; a terminal 
voltage detecting section capable of measuring a 
voltage across terminals of the secondary cell; a 

10 parameter estimating section that calculates an 

adaptive digital filtering using a cell model in a 
continuous time series shown in an equation (1) and 
estimates all e-f — parameters at one time, the 
parameters corresponding to an open-circuit voltage 

15 Vq, which is an offset term of the equation_^ (1) and 
coefficients of A(s), B(s), and C(s)^ which are 
transient terms; and a charge rate estimating section 
that estimates the charge rate from a previously 
derived relationship between a previously — derived an 

20 open-circuit voltage and a charge rate of the 

secondary cell and the open-circuit voltage Vo and 

th-e — charge — rate — &0€ — using — trh-e — open — circuit — voltage — V-e-, 

V = -^^•/+ — ^ — •Fn (l)f wherein s denotes a 

A(s) C(s) ' 

25 Laplace transform operator, A(s), B(s), and C(s) 
denote poly-nominal functions of s. 

[0006] According to another aspect of the present 

invention, there is provided a charge rate 
estimating method for a secondary cell, comprising: 
30 measuring a current flowing through the secondary 
cell; measuring a voltage across terminals of the 
secondary cell; calculating an adaptive digital 
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filtering using a cell model in a continuous time 
series shown in an equation (1); estimating all 
parameters at one time, the parameters corresponding 
to an open-circuit voltage Vp, which is an offset 
5 term of the equation (1)^ and coefficients of A(s), 
B(s), and C ( s ) _^ which are transient terms; and 
estimating the charge rate from a previously derived 
relationship between an open-circuit voltage and a 
charge rate of secondary cell^ and the e — proviouoly 
10 derived — open — circuit — voltage — V^ ^ and — t^i-e — charge — rate 
SQC — using — t4^-e — open-circuit voltage Vq, 

V = — ^ — •Vn (1)/ wherein s denotes a 

A(s) C(s) ° 

Laplace transform operator, A(s), B(s), and C(s) 
denote poly-nominal functions of s, 

15 [0007] According to a still another object of the 

present invention, there is provided a charge rate 
estimating method for a secondary cell, comprising: 
measuring a current I(k) flowing through the 
secondary cell; measuring a terminal voltage V(k) 

20 across the secondary cell; storing the terminal 

voltage V(k) when a current is zeroed as an initial 
value of the terminal voltage AV ( k ) = V{k) - V_ini; 
determining instantaneous current values Io(k), Ii(k), 
and l3(k) and instantaneous terminal voltages Vi(k), 

25 V2(k), and V3(k) from an equation (19), 

10 =—!—•/, 

11 =-^#/, Vi = •F, 

12 = •// V2 = •V r 
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10 



I3 = 



V3 = 



and 



1 



1 



C?,(5) 

(19), wherein pi denotes a 



G,(^) {pi»s+iy 

constant determining a responsive characteristic of 
Gi(s); substituting the instantaneous current values 
Io(k), Ii(k), 1 2 ( k ) ^ a nd l3(k) and the instantaneous 
terminal voltages Vi(k), V2(k), and V3(k) into an 
equa t ion (18), 



Y(k) 
e (k) 

P (k) 

>-i(k) = 



l + Zj(k)»oi^(k)» P(Jc - 1) • coik) 

e(k -1) - Y(k) (k-1) •CO (k) • [co"^ (k) •e (k -1> -y(k)] 



1 



\P{k-\) 



X,{k) 
fy^trace{p\k) ] 



A3 {k) • P(k - 1) • (o{k) •co'^ik)* P(k - 1) 
\-^X^{k)»co'^{k)*P(k-\)»co(k) 



P\k) 



i- 



trace 



ru 



^ ^1 ' 

trace 



trace 



Ip'm } 



U 



Yu 
{P'(k) } 



<^i< 



trace 



[P\k) ] 



Yl 



trace 



Yl 



Yl 



< u 



15 



20 



(18 ) , 

wherein 9(k) denotes a parameter estimated value at a 

time point of k (k = 0, 1, 2, 3 ), Xx, XaCk), yu , 

and yL denote initial set value, b < Xi < 1, 0 < A,3(k) 
< 00. P(0) is a sufficiently large value, 9(0) 
provides an initial value which is non-zero but very 
sufficiently small value, trace{P} means a trace of 
matrix P, wherein y(k) = Vi(k) 
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co^(k) = [V3(k) V2(Jc) l3(k) l2(k) Ii(k) 

lo (k) ] 
e (k) = 



substituting a, c, e, and f in the parameter 

5 estimated value G(k) into and equation (22) to 
calculate Vq * which is an alternate of Vo which 
corresponds to a variation AVq ( k ) of the open-circuit 
voltage estimated value from a time at which the 
estimated calculation start is carried out; 

10 Vo = U ^ X • ^0 = a • Vg + b • Vs +V4 -c • le - d • I5 - 

G2 W 

e • I4 (22); and calculating an open-circuit 

voltage estimated value Vo(k) according the variation 
AVo(k) of the open-circuit voltage estimated value 
and the terminal voltage initial value V_ini. 
15 [0008] This summary of the invention does not 

. necessarily describe all necessary features so that 
the invention may also be a sub-combination of these 
described features. 

20 BRIEF DESCRIPTION OF THE DRAWINGS: 

[0009] Fig. 1 is a functional block diagram of an 

apparatus for estimating a charge rate (SOC) of a 
secondary cell in a preferred embodiment according to 
the present invention. 

25 [0010] Fig. 2 is a specific circuit block diagram 

of the apparatus for estimating the charge rate ( SOC ) 
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of the secondary cell in the preferred embodiment 
according to the present invention. 
[0011] Fig. 3 is a model view representing an 
equivalent circuit model of the secondary cell . 
5 [0012] Fig. 4 is a correlation map representing a 
correlation between an open-circuit voltage and a 
charge rate (SOC) . 

[0013] Fig. 5 is an operational flowchart for 

explaining an operation of a microcomputer of a 
10 battery controller of the charge rate estimating 

apparatus in the first preferred embodiment shown in 
Fig . 1 . 

[0014] Figs. 6A, 6B, 6C, 6D, 6E, 6F^ 6G , 6H and 61 

are characteristic graphs representing results of 
15 simulations of current, voltages, and various 

parameters in a case of the charge rate estimating 
apparatus in the embodiment shown in Fig. 1 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS: 
[0015] Reference will hereinafter be made to the 
drawings in order to facilitate a better 
understanding of the present invention. 
[0016] Fig. 1 shows a functional block diagram of 

charge rate estimating apparatus in a first preferred 
embodiment according to the present invention. In Fig. 
1, a reference numeral 1 denotes a parameter 
estimating section based on a cell model with an 
open-circuit voltage Vq (k) as an offset term. In 
addition, a reference numeral 2 denotes a open- 
circuit voltage calculating section to calculate 
open-circuit voltage Vq (k), and a reference numeral 
3 denotes a charge rate estimating section that 
calculate the charge rate from the open-circuit 
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voltage. In addition, a reference numeral 4 denotes 
a current I measuring block to detect current I(k) 
which is charged and discharged into and from the 
cell, and a reference numeral 5 denotes s terminal 
5 voltage of the cell to measure the terminal voltage 
V(k) . 

[0017] Fig. 2 shows a block diagram representing a 

specific structure of the charge rate estimating 
apparatus in the first embodiment. In this embodiment, 

10 a load such as a motor is driven with the secondary 
cell and the charge rate estimating apparatus is 
mounted in a system to charge the secondary cell with 
a regenerative power of the motor (load) . In Fig. 2, 
a reference numeral 10 denotes a secondary cell 

15 (simply called, a cell), a reference numeral 20 
denotes a load such as a DC motor, a reference 
numeral 30 denotes a battery controller (electronic 
control unit) to estimate the charge rate (charge 
state) of the cell having a microcomputer including a 

20 ROM (Read Only Memory) , a RAM (Random Access Memory) , 
a CPU (central Processing Unit), and Input/Output 
Interface and other electronic circuits. A reference 
numeral 40 denotes a current meter to detect a 
current which is charged into or discharged from the 

25 cell, a reference numeral 50 denotes a voltage meter 
to detect the terminal voltage of the cell, a 
reference numeral 60 denotes a temperature meter to 
detect a temperature of the cell. These meters are 
connected to battery controller 30. Battery 

30 controller 30 corresponds to parts of parameter 

estimating section 1, an open-circuit voltage Vq (k) 
and a charge rate estimating section 3. Current 
meter 40 corresponds to current I(k) measuring 
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section and voltage meter 50 correspond to terminal 
voltage V(k) measuring section 5. 

[0018] First, a " cell model " used in the first 

embodiment will be described below. Fig. 3 is an 
5 equivalent circuit representing an equivalent circuit 
model of the secondary cell. The equivalent circuit 
model of the secondary cell can be represented by the 
following equation (7) (= equation (6)). 



10 V = ^ — r-^*^ +7;^ 7 Vo (7) . 

7; • 5- + 1 73 • 5 H- 1 

In equation (7), a model input is a 
current I [A] (a positive value represents a charge 
and a negative value represents a discharge), a model 
output is a terminal voltage V [V] , an open-circuit 

15 voltage is Vo, K denotes an internal resistance, Ti 

through T3 denote time constants (Ti 5t T2 ^ T3, Ti« T3) 
and s denotes a Laplace transform operator. 
[0019] In this model based on equation of (7) is a 
reduction model (first order) in which a positive 

20 pole and a negative pole are not specially separated 
from each other. However, it is possible to represent 
a charge-discharge characteristic of an actual cell 
relatively easily. Equation (7), in equation. (1) of V 
= B(s)/A(s)*I + l/C(s)»Vo (1), A(s) = Ti • s + 1, 

25 B(s) = K • (T2-)- • s + !_)_, C(s) = T3 • s + 1 . 

[0020] Hereinafter, a deviation from the cell model 
based on equation (7) to an adaptive digital filter 
will first be described below. Open-circuit voltage 
Vo can be described by an equation (8), supposing 

30 that a value of a current I multiplied with a 
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variable efficiency of A is integrated from a certain 
initial state . 

A 

That is to say, Vq = — •/ (8). 

s 

It is noted that equation (8) corresponds to 
5 a replacement of h recited in equation (2), viz., Vq 
= h/s • I with efficiency of A. 

[0021] If equation (8) is substituted into equation 
(7), equation (9) is resulted. 



10 



15 



20 



25 



Vo = 



1 



-•/ (9) 

s 



t; • 5 + 1 

Equation (9) corresponds to equation (3) (V 




B(s)^ 1 



[^A(s) C(s) s 



- .1 = 



S^Bis)^C(s)^i'h^Ais) 
S^A(s)^Cis) 



• I 



'(3)). For A{s), B(s), and C(s) in equation (3), 

the following equations are substituted into equation 

(9) in the same way as the case of equation (7) . 

A ( s ) = Ti • s + 1 , 
B (s) = K • {T2*s + 1) 

C(s) = Ta's + 1. In other words, equation 
(3) is a generalized equation and this application to 
a first order model is equation (9). If equation (9) 
is arranged, an equation of (10) is given. 
S •(Ti*s + 1) (T3»s + l)»V = KMT2»s + 1) {T3*s + l)s* I + 
A^(Ti*s + 1)*I 

{Tl^ T3^ + (Ti + T3)»S^ + s} • V = {K*T2-T3#S^ + K • 

{T2 + T3)'s^ + (K + A • Ti)*s + A} • I 

•^ + e*s + f)»I 

(10) . It is noted that, in the last equation of 
equation (10), parameters are rewritten as follows: 



(a»s^ + b»s^ + s)»V = (c*s^ + d^s^ 
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a = Ti»T3/ b = Ti + T3/ c = K • T2 • T3, d = K • {^^T_2 + 

T3) , e == K + A • Ti, and f = A (11). 

If a stable low pass filter Gi(s) is introduced into 
both sides of equation (10) and arranged, the 
following equation (12) is given. 



(a-s-^ H- b^s"" + s ) •V = (c^s-" + d»s^ + e«s + f)»I 



(12). In details, in equation (10), on the 

contrary of equation (7), if Ti • s + 1 = A(s), 

(T2*s + 1) = B(s), and T3*s + 1 = C{s) are substituted 

10 into equation (10), this is given as: s«A ( s ) •C ( s ) •V = 
B is) •C is) •s* I + A»A(s)»I. This is rearranged as 
follows: s^A{s) •C (s) •V = [B(s)*C(s)»s [[-I]] + 

A»A(s)]*I (12)'. If, the low pass filter (LPF), 

Gi{s) is introduced into both sides of equation (12) 

15 an equation (4) is given. 

That is to say, — ^^•F = — — ^#7 

(4) . It is noted that s denotes the Laplace 
transform operator, A(s), B(s), and C(s) denote a 

20 poly-nominal function of s, h denotes a variable, and 
1/Gi(s) denotes a transfer function having a low pass 
filter characteristic. That is to say, equation (4) 
is the generalized function, equation (12) is the 
application of equation (4) to the first order model. 

25 [0022] Current I and terminal voltage V which can 

actually be measured are processed by means of a low 
pass filter (LPF) and a band pass filter (EPF) are 
defined in the following equations (13), provided 
that pi denotes a constant to determine a responsive 
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characteristic of Gl(s) and is determined according 
to a designer's desire. 

1 

10 = •/ 

G,(s) 

11 = —^•I, Vi = —^•V, 

12 = — — V2 = — •¥ , 

G,{s) G,(5) 

13 = — — •/, V3 = — •¥ , 

G,{s) G,is) 



1 1 



(13) 



If equation (12) is rewritten using the variables 
10 shown in equations (13), equations (14) are 

represented and, if deformed, the following equation 
(15) is given . 

a'Vs + b«V2 + Vi = c«l3 + d»l2 + e»Ii + 

f • lo 

15 Vi = -a«V3 - b»V2 + c»l3 + d»l2 + e • Ii + 

f»Io (14) . 



20 



Vi = [V3 V2 I3 



lo] • [ [ = ] ] 



(15) . 



— a 
-b 
c 
d 

e 

J . 

Equation (15) is a product-sum equation of measurable 
values and unknown parameters. Hence, a standard 
(general) type (equation (16)) of the adaptive 
digital filter is coincident with equation (15). 
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It is noted that co"^ means a transposed vector in 
which a row and column of a vector co are mutually 
exchanged . 

y = ©'^•G (16) - It is noted that y, 



CO 



and 9 can be expressed in the following equation (17) 
in equation (16) described above. 



Y = Vi, (0^ = [V: 



V2 



lo] . e = 



(17) 



- a 
-b 
c 
d 
e 

f . 

Hence, if a signal filter processed for current I and 
terminal voltage V is used in a digital filter 
10 process calculation, unknown parameter vector 9 can 
be estimated . 

[0023] In this embodiment, " a both-limitation 

trace gain method is used which improves a logical 
demerit of a simple " an adaptive digital filter by 

15 means of a least square method " such that once the 
estimated value is converged, an accurate estimation 
cannot be made any more even if the parameters are 
changed. A parameter estimating algorithm to estimate 
unknown parameter vector 9 with equation (16) as a 

20 prerequisite is as shown in an equation (18) . It is 
noted that the parameter estimated value at a time 
point of k is 9(k) . 



25 



Y(k) = 



l-hA,(k)^(D^ (k) • P(k - 1) • a)(k) 



9(k) = 9(k -1) - Y ( k) ( k-1 ) ( k) • [CD^ ( k) •e ( k -1> -y(k)] 
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P (k) = 



1 



Pik-l)- 



tra ce\p'ik) } ^ 



A3 (k) • P(k - 1) * o)(k) • a/ {K) • Pjk - 1) 
\-k-X^(k)»a)''ik)»P{k-\)»co{k) 

trace{p\k) ] 



It 



trace {P (k) } _ trace 



Yu 

trace{p'ik) } 
Yl 



■<A,i< 



[Pik-) } 



trace 



Yl 



Yl 



and Yl 
X3 {k) < 00, 



(18) . 

In equations (18), Xi, Xsik), jui 
denote initial set value, h < Xi < 1 , 0 < 
P{0) is a sufficiently large value, 9(0) provides an 

10 initial value which is non-zero but very sufficiently 
small value. In addition, trace{P} means a trace of 
matrix P. As described above, the derivation of the 
adaptive digital filter from cell model. 
[0024] Fig. 5 shows an operational flowchart 

15 carrying out the microcomputer of battery controller 
30. A routine shown in 5 is carried out for each 
constant period of time To . For example, I(k) is the 
present value and I(k -1) means a one previous value 
of I(k) . At a step SIO, battery controller 30 

20 measures current I{k) and I(k -1) means one previous 
value of I{k). At -s—a step S20, battery controller 30 
carries out a turn on-and-off determination of an 
interrupt relay of the secondary cell. That is to 
say, battery controller 30 performs the on-and-off 

25 control of the interrupt relay of the secondary cell. 
When a relay is turned off (current 1=0), the 
routine goes to a step S30. During the engagement of 
the relay, the routine goes to a step S40. At step 
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S30, when the relay is engaged, the routine goes to a 
step S540. At step S530, battery controller 30 
serves to store terminal voltage V(k) to as an 
initial value of the terminal voltage V_ini. At a 
5 step S40, battery controller 30 calculates a 

differential value AV ( k ) of the terminal voltage. 
AV(k) = V(k) - V__ini . This is because the initial 
value of the estimation parameter in the adaptive 
digital filter is 0 so that the estimation parameter 

10 does not converge during the estimation calculation 
start time. Thus, all of inputs are zeroed. During 
the input being all zeroed. During the relay 
interruption, step S30 have been passed and the 
estimation parameters are remains initial state since 

15 1=0 and the estimation parameter remains alive. 

[0025] At step S50, a low pass filtering or band 

pass filtering are carried out the current I(k) and 
terminal voltage difference value AV ( k ) on the basis 
of equation (13) - Io(k) through l3(k) and Vi(k) 

20 through V3(k) are calculated from equation (19) . In 

this case, in order to improve an estimation accuracy 
of the parameter estimation algorithm of equation 
(18), a responsive characteristic of low pass filter 
Gi(s) is set to be slow so as to reduce observation 

25 noises. However, if the characteristic is quicker 

than a response characteristic of the secondary cell 
(a rough value of time constant Ti is known), each 
parameter of the electric cell model cannot 
accurately be estimated. It is noted that pi recited 

30 in equation (19) denotes a constant determined 

according to the responsive characteristic of Gi(s) . 
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15 



[0026] 



1 

s 

„2 



10 = 

11 = 

12 = 

13 = 
1 ^ 1 



G,(5) 



Vi = 



V2 = 



V3 = 



G,(s) 



V , 



G,{s) 



• V , and 



(19) 



At a step S60, Io(k) through l3(k) calculated at step 
S50 and Vi(k) through V3(k) are substituted into 
equation (18) . Then, the parameter estimation 
10 algorithm in the adaptive digital filter, viz., 
equation (18) is executed to calculate parameter 
estimated value 9(k) . y(k) , to''(k) , and e(k) are shown 
in equation (20). 

y ( k) = Vi ( k) 

co^(k) = [V3(k) V2(k) l3(k) l2(k) Ii(k) 

lo (k) ] 



e (k) = 



-bik) 

cdk) 

d(k) 

e(k) 

f(k) 



(20) 



At a step S70, a through e of parameter estimated 
20 value 9(k) calculated at step S60 are substituted 

into the following equation (22) in which the above- 
described cell model equation (7) is deformed to 
calculate Vo ' which is an alternative to open-circuit 
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voltage Vo. Since the variation in open-circuit 
voltage Vo is smooth, Vq ' can be used alternatively- 
It is noted that the derivation herein is a variation 
AVo(k) of the open-circuit voltage from the estimated 
5 calculation start time. 

[0027] It is noted that an equation of [l/Cl(s)]I 

in equation (21) is replaced with an equation (24)' 
corresponds to equation (22) . It is also noted that, 
in the derivation of equation (22), K in equation 
10 (21) is strictly different from e in equation (21). 
However, since, physically, K » A • Ti, e is 
approximated to K (e = K) . Then, each coefficient a 
through e in equation (22) is the contents shown in 
equation (23). 

15 •Fo = V - ^ 

(Ti«s + !)• Vo =(Ti#s + 1) (T3«s H- 1)V - 
K» (T2*s + 1 ) (T3»s + 1 ) #1 

(Ti#s + !)• Vo = {Ti*T3*s^ + (Ti + T3)»s + 

1}»V 

20 - {K»T2*T3*S^ + K«(T2 +T3)*S + K}»I 

(21) , 

V'o = * ^) ^ = ^•Ve + b^Vs +V4 -c^Ie - d^Is - e« I4 - 

Gi is) 



25 



(22) . 

[0028] It is noted that a = Ti • T3, b =, Ti + T3, c 

= K •(T2 + T3}., d = K •(T2 + T3) , e = K + A«Ti = K-- 

- (23) . 
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10 



15 



20 



25 



[0029] 



l4 = 



Is = 



V4 = 



1 



1 



Vs = 



•V , 



G^{s) 

l6 = 



P2»s+\ r, '-J + 1 ' 



• / , and Ve = 



(24) 



[0030] p2 recited in equations (24) denote a 

constant to determine a responsive characteristic of 
G2 ( s ) . Ti of the ceil parameter is known to be 
several seconds. Hence, T'l in equation (24) is set 
to be approximated value to Ti. Thereby, since (Ti» s 
+ 1 ) which remains in a numerator of equation (22) 
can be compensated, the estimation accuracy of open- 
circuit voltage Vq can be improved. It is noted that 
equation (21) corresponds to equation (5) , That is to 
say, equation (21) is derived from (Ti • s +!)• Vo = 
(Ti • s + 1) (T3 • s + 1)»V - K«(T2« s + 1) (T3» s +!)• 
(T3 • s + !)• I. If the following three equations are 
substituted into the above-described deformation of 



equation (21). 



s+l=A(s), K» (T2 



s + 1 ) 



B(s), and T3 • s + 1 = C(s). That is to say, A(s)« Vo 
= A(s)« C(s) • V - B(s)» C(s)» I. If this is 
rearranged, this results in Vq = C(s) • v - B(s) • 
C(s) • I/A(s), Vo = C(s)»[V - B(s) • I/A(s)] If low 
pass filter G2(s) is introduced into both sides of 
this equation, this results in equation (5) . In 
details, equation (5) is a generalization equation 
and the application of equation (5) to the first 
order model is equation (2) . 
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[0031] At a step S80, battery controller 30 adds 

the open-circuit voltage initial value, i.e., 
terminal voltage initial value V__ini to a variation 
AVo(k) of open-circuit voltage Vq so as to obtain 
5 open-circuit voltage estimated value Vo(k:) from the 
following equation (25) . 

Vo (k) = A Vo (k) + V_ini (25). 

[0032] At a step S90, battery controller 30 

calculates the charge rate SOC(k) from open-circuit 

10 voltage Vq (k) calculated at step S80 using a 

correlation map of the open-circuit voltage versus 
the charge rate as shown in Fig. 4. It is noted that, 
in Fig. 4, Vl denotes the open-circuit voltage 
corresponding to SOC = 0 % and Vh denotes the open- 

15 circuit voltage corresponding to SOC = 100 %. At a 

step SlOO, battery controller 30 stores the necessary 
numerical values needed in the subsequent calculation 
and the present routine is ended. As described above, 
an operation of the apparatus for estimating the 

20 charge rate of the secondary cell has been described. 
[0033] (1) As described above, a relationship 

from among current I of the secondary cell and 
terminal voltage V thereof, and the open-circuit 
voltage Vq is structured in transfer function that as 

25 in the general equation (1), that in the preferred 
embodiment, equation (7) (= equation (6). Hence, it 
is made possible to apply an adaptive digital filter 
such as a least square method (well known estimation 
algorithm) . Consequently, it becomes possible to 

30 estimate parameters in equations (viz., open-circuit 
voltage Vq which is an offset term and poly-nominal 
equations A(s), B(s), and C(s)) in a form of a batch 
processing. These parameters are largely affected by 
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the charge rate, a surrounding temperature, and a 
deterioration and varied instantaneously. It is 
possible to sequentially estimate the adaptive 
digital filter with good accuracy. Then, if a unique 
5 correlation between the open-circuit voltage Vq and 
the charge rate as shown in Fig. 4 are stored, the 
estimated open-circuit voltage can be converted to 
the charge rate. Hence, it is possible to 
sequentially estimate the charge rate in the same way 

10 as the parameters described above. 

[0034] (2) In a case where the equation (1) which 

is the relationship equation of current I and 
terminal voltage V of the secondary cell is 
approximated to equation (4), the equation such that 

15 no offset term is included (viz., the open-circuit 

voltage Vq), a product-and-addi tion equation between 
a measurable current I which is filter processed and 
a terminal voltage V which is filter processed and 
unknown parameter (coefficient parameters of poly- 

20 nominal equations A(s), B(s), and C(s) and h) is 
obtained. A normally available adaptive digital 
filter (the least mean square method and well known 
parameter estimation algorithm) can directly be 
applied in a continuous time series. 

25 [0035] As a result of this, the unknown parameters 

can be estimated in the batch processing manner and 
the estimated parameter h is substituted into 
equation (2), the estimated value of open-circuit 
voltage Vq can easily be calculated. All of these 

30 parameters are varied instantaneously, the adaptive 
digital filter can serve to estimate the charge rate 
at any time with a high accuracy. Since a constant 
relationship between open-circuit voltage Vq and the 
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charge rate J_SOCJ^ is established as shown in Fig, A, 
if this relationship is previously stored, the charge 
rate {JSOC)__ can be estimated from the estimated value 
of open-circuit voltage Vq • 
5 [0036] Figs. 6A through 61 integrally shows 

signal timing charts with current I and terminal 
voltage V inputted into adaptive digital filter and 
representing results of simulation graphs when each 
parameter is estimated. As far as a time constant of 

10 a first order delay in equation (6) is concerned, Ti 
< To. Since all parameters a through f (refer to 
equation (11)) are favorably estimated, the estimated 
value of open-circuit voltage Vq can be said to be 
well coincident with a real value. 

15 [0037] It is noted that, in Fig. 6C which 

indicates the open-circuit voltage, a reason that a 
right side second term of equation (6) is described 
is to indicate that the open-circuit voltage 
estimated value is coincident with a real value 

20 almost without delay in spite of the fact that a 
late term of time constant T3 is measured on the 
terminal voltage inputted into the adaptive filter. 
In details, since the parameter estimation with the 
cell model formatted adaptive digital filter in 

25 equation (6), all of parameters a through f can 

favorably be estimated and the estimated value of 
open-circuit voltage Vq is well coincident with a 
real value . 

[0038] (3) In addition, as described in item (2), 

30 in the structure in which the open-circuit voltage Vq 
is calculated from equation (2), the integration 
occurs before a value at which estimated value h is 
converged to the real value, its error cannot be 
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eliminated. However, in the structure in which 
equation (5) in which the integration is not included, 
the error before the parameter estimated value is 
converged into the real value does not give an 
5 influence after the convergence. 

[0039] It will be appreciated that, in part of ® 

in Fig. 61, before estimated value f is converged 
into a real value, an erroneous estimation is carried 
out only at momentarily. In equation (2), this value 

10 is also integrated so that the error is not 

eliminated. However, the error is not eliminated 
since even this value is integrated. However, in the 
structure using equation (5), open-circuit voltage Vq 
is calculated from the equation in which the 

15 integration is not included. Hence, after the 

parameter estimated value is converged into the real 
time, this erroneous estimation portion is eliminated. 
[0040] (4) Furthermore, in a case where equation 
(6) is used in place of equation (1), a calculation 

20 time and program capacity can be suppressed to a 

minimum while having the above-described advantages. 
[0040] [0041] The entire contents of a Japanese 

Patent Application No. 2002-340803 (filed in Japan on 
November 25, 2002) are herein incorporated by 

25 reference. The scope of the invention is defined with 
reference to the following claims. 
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